Cosmic Flows on 100 h Mpc scales 
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To study galactic motions on the largest available scales, we require bulk flow moments whose 
window functions have as narrow a peak as possible and having as small an amplitude as 
possible outside the peak. Typically the moments found using the maximum likelihood esti- 
mate weights do not meet these criteria. We present a new method for calculating weights 
for moments that essentially allow us to "design" the moment's window function, subject, of 
^"v course, to the distribution and uncertainties of the available data. 
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On scales that are small compared to the Hubble radius, galaxy motions are manifest in 
deviations from the idealized isotropic cosmological expansion 



^ cz = H Q r + f • [v(r) - v(0)] . (1) 

Redshift-distance samples, obtained from peculiar velocity surveys, allow us to determine the 
radial (line-of-sight) component of the peculiar velocity of each galaxy: 

v(r) = r • v(r) = cz — Hqt . (2) 



Galaxies trace the large-scale linear velocity field v(r) which is described by a Gaussian 
random field that is completely defined, in Fourier space, by its velocity power spectrum P v (k). 
£\\ [ The Fourier Transform of the line-of-sight velocity is 
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r • v (r) = — -„ 
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which provides us with the velocity power spectrum 



r • v (r) = / d 3 kf-k ^(k)e ik ' r , (3) 
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where <5d is the Dirac delta function. In linear theory, the velocity power spectrum is related to 



< v(k)v*(k') >= (27rfP v (k)S D (k - k') , (4) 
f 

the density power spectrum 



Pv(k) = ^rf 2 (n mfl n A )P(k) , (5) 

where H is the Hubble constant and / is the rate of perturbation growth at the present. O m and 
£l\ are the matter and dark energy density parameters today. The power spectrum provides a 
complete statistical description of the linear peculiar velocity field. 

For our power spectrum model, we follow Eisenstein and Hu (1998) in writing 

P{k) cx a 8 n° m 6 k n T 2 {k/T) , (6) 

where T is the transfer function, T is the power spectrum "shape" parameter and as usual, 
we normalize the power spectrum using the parameter erg, the amplitude of matter density 
perturbations on the scale of 8 /i~ 1 Mpc. 

To find the maximum likelihood likelihood parameters from peculiar velocity surveys we 
start with catalog of peculiar velocities galaxies, labeled by an index n, positions r n , estimates 



of the line-of-sight peculiar velocities S n and uncertainties a n . We assume that the observa- 
tional errors are Gaussian distributed. We then model the velocity field as a uniform streaming 
motion, or bulk flow (BF), denoted by U, about which there are random motions drawn from 
a Gaussian distribution with a 1-D velocity dispersion cr*. Then the likelihood function for the 
BF components is 

where i go from 1 to 3 to specify the BF components. While this has been done previously for 
maximum likelihood estimates (MLE) of the bulk flow moments, the advantage of these new 
moments is that they have been designed to be sensitive only to scales of order 100ft. _1 Mpc. 
Thus we will be able to probe these scales without having to worry about the influence of 
smaller scales. Further, by isolating the very large scale motions we will be able to put stronger 
constraints on power spectrum parameters. 

Maximum likelihood solution for BF moments is 



where 



The measured peculiar velocity of galaxy n 

S n = r n ,iVi(r n ) +e n , (10) 

where e n is a Gaussian with zero mean and variance o~\ + o~1 . The theoretical covariance matrix 
for the BF components 

R ab = R^+R$. (11) 
The first term is given as an integral over the matter fluctuation power spectrum, P(k), 



Kab ~ 2tt 2 



/ dk W 2 ab (k)P(k), (12) 
Jo 



where the angle-averaged tensor window function is 

Wab(k) = w a,nW b , m f ^ (? n ■ k f m ■ k) exp (ik • (r n - r m )) . (13) 

n,m 

In order to do better than the MLE solutions, we redesign the window functions by requiring 
that they have narrow peaks and small amplitude outside the peak. We start with the BF MLE 
weights 

n u n ~T u * 



which depends on the spatial distribution and errors. Now consider an ideal survey with very 
large number of points, isotropic distribution and a Gaussian falloff n(r) cx exp(— r 2 /2Rj) where 
Ri is the depth of the survey. 

The moments are specified by the weights 



Ui — Wi^ n S n 

n 



(15) 
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Fig. 1: The window functions of the bulk flow com- 
ponent for Rf = 50?t _1 Mpc for the all the catalogs we 
considered. The thick (thin) lines are the window func- 
tions for the MV (MLE) bulk flow components, the 
x,y,z-component are dash-dot, short dash, long dash 
lines respectively. The Thick solid line is the ideal win- 
dow function. 




0.2 0.4 0.6 0.80 0.2 0.4 0.6 0.6 



0.2 0.4 0.6 0.6 



Fig. 2: The T likelihood function for the composite 
catalogs. 



that minimize the variance < (v,i — Ui) 2 >. We expand the variance 

< (Ui - Ui) 2 >= W i)n Wi, m < S n S m > + < U 2 > -2 Y w i,n < UiS n > , (16) 



since the measurement error included in S n is uncorrelated with the bulk flow Ui. For BF 
moments we can impose the constraint that 
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(17) 



We now minimize this expression with respect to Wi^ n subject to the constraint which we 
enforce using a Lagrange multiplier to get 



Wj 



(S n Ui, 



(18) 



We solve this to get the minimum variance (MV) weights. 

To apply this formalism, we used both cluster and galaxy peculiar velocity surveys. Specif- 
ically, we used the Willick cluster survey^, the SClE an d 802© (Tully-Fisher clusters); The 
EFAR survey^ (a Fundamental Pla ne clu ster survey); the SN are Type la supernovae from t he 
compilation of Tonry et al. 

B S MAC™ 

(a Fundamental Plane cluster survey); SFI+-P^^ (a 
Tully-Fisher based survey of ~ 2700 field galaxies and ~ 730 groups). 

The minimum variance, or MV, weights were calculated for the bulk flow component mo- 
ments using the method described above for each of our catalogs, as well as for three composite 
catalogs, SFI++, consisting of both fields and groups from the catalog, OTHER, consisting of 
all the surveys except for the SFI++ field galaxy and groups catalogs, and Composite, consisting 
of all of the catalogs combined. The OTHER catalog was included due to the fact that, given 
their large size, the SFI++ catalogs tend to dominate any composite catalog in which they are 
included. Here we will show results from a deep survey, Rj = 50 /i _1 Mpc. 
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Fig. 3: The bulk flow of the large catalogs. The triangles (points) are the BF moments with 
Rj = 20/i _1 Mpc = 50/i _1 Mpc ). There is a consistent and robust flow exhibited in all 

catalogs to the negative Galactic y direction (upper right panel) for the large scale analysis which 
is reflected in the BF magnitude (lower right panel) 



In Figure 1 we show the window functions of the MV bulk flow component moments. For 
comparison, we also include the ideal window functions as well as those for the MLE moments 
for each survey. As expected, the match between the window functions for the MV moments 
and the ideal is best for the large surveys and those with small measurement error and similar 
distribution to the ideal survey. For the sparse, noisy surveys, the window functions for the 
MV moments are not very different than those of the MLE moments, differing mostly in the 
amplitude of the tail of the window function for large k. 

In Figure 2. we show the likelihood function of the "shape" parameter T for the large 
catalogs. We see clearly that T is smaller for large Ri than for smaller Rj. In Figure 3. we 
show the BF velocity for the composite surveys. We see that for the shallow window functions, 
we get small flows, ~ 150 km/s. On large scales we get very large flows ~ 400 km/s. This 
unexpected flow which is exhibited in all the catalogs we analyzed, require a very steep power 
spectrum which leads to small T (see figure 2.) 
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